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Abstract 

The data on the total cross section atot{G'^G~ ^e~^e~bb) measured at LEP2 represent a serious challenge 
for perturbative QCD. In order to understand the origin of the discrepancy between the data and QCD 
I calculations, we investigate the dependence of four contributions to this cross section on 77 collision energy. 
CS| ■ As the reliability of the existing calculations of atot{e'^e~ — >e"'"e~66) depends, among other things, on the 
^ I stability of the calculations of the cross section atot (77 bb) with respect to variations of the renormalization 
O ■ and factorization scales, we investigate this aspect in detail. We show that in most of the region relevant 
for the LEP2 data the existing QCD calculations of the cross section atoti'JJ — > bb) do not exhibit a region 
■ of local stability. Possible source of this instability is suggested and its phenomenological implications for 
I understanding the LEP2 data are discussed. 

(N , 

_^ '. 1 Introduction 

'NT 

I> 

' Heavy quark production in hard collisions of hadrons, leptons and photons has been considered as a clean 
test of perturbative QCD. It has therefore come as a surprise that the first data on the bb production in pp 
collisions at the Tevatron 1 7P collisions at HERA 3 , 4 and 77 collisions at LEP2 (^BHI have turned out 
to lie significantly and systematically above theoretical calculations. The disagreement between data pHO] 
and theory |71|Hlini was particularly puzzling for the collisions of two quasireal photons at LEP2. 

' ^ • The arrival of new data on bb production in cp collisions at HERA ^U], shown in the left part of Fig. ^as 
I I solid squares, have further complicated the situation. In the range of moderate ~ 80 GeV^ the new ZEUS 
data are in reasonable agreement with NLO QCD predictions and also in the photoproduction region the 
(— I I excess of the new data over theory is substantially smaller then that in the older data. As a result, there is 
now an inconsistency between new ZEUS and older HI results 131 for moderate Q^, but the situation remains 
unclear also in the photoproduction region. For pp collisions the progress on the theoretical side (llil2) has 
significantly reduced the discrepancy observed at the Tevatron. 

^ . On the other hand, the problem of understanding the bb production in 77 collisions remains. The preliminary 

DELPHI data presented this spring at PHOTON 2003 conference JHI and reproduced in Fig. ^ are in striking 
agreement with the older L3 and OPAL data. The central values of all three experiments are almost identical 
which strongly supports the reliability of these measurements. Contrary to the case of analogous discrepancy 
in antiproton-proton collisions at the Tevatron, there have been few theoretical suggestions how to explain the 
sizable excess of data over current theory in 77 collisions. Neither the use of unintegrated parton distribution 
functions jl4| . nor the production of supersymmetric particles '15', proposed for explaining an analogous 
excess in antiproton-proton collisions, are of much help for LEP2 data, primarily because of low 77 energies 
involved. Quite recently, however, this discrepancy has been interpreted as an evidence for integer quark 
charges We will come back to this suggestion in Section 4.1. 

In ^21 we have investigated the sensitivity of QCD calculations of (Jtotipp bb; S, M, ^) to the variation of 
the renormalization and factorization scales fi and M. In particular we have argued that in order to arrive 
at locally stable results JJj these two scales must be kept independent. We have furthermore shown that 
in the Tevatron energy range the position of the saddle point of the cross section atot {pp — > bb; S, M, yit) 
lies far away from the "diagonal" ji = M used in all existing calculations. Using the NLO prediction at 
the saddle point instead of the conventional choice ^ — M — mi, enhances the theoretical prediction in the 
Tevatron energy range by a factor of about 2, which may help explaining the excess of data over NLO QCD 
predictions. 

In this paper similar analysis is performed for 77 collisions at the total centre of mass energies W relevant 
for existing LEP2 data. The specific features of the theoretical description of QQ production in 77 collisions 
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Figure 1: The current situation with data on hh production in ep and e+e colhsions, including the most 
recent data of ZEUS [101 and DELPHI [H]. 



have been discussed in JH|- However, as all three experiments at LEP2 have measured merely an integral 
over the cross section (7(77 bb, W) weighted by the product of photon fluxes inside the beam electrons 
and positrons, it is important to understand the W^-dependence of the four individual contributions to it. 

The paper is organized as follows. The basic facts and formulae relevant for the quantitative investigation 
of renormalization and factorization scale dependence of finite order QCD approximations are collected in 
Section 2. This is followed in Section 3 by the discussion of the general form of dtot (77 bb; W, M, /i). In 
Section 4 the I^-dependence of the four contributions to the cross section (Jtot{e'^e~ — >e~'"e~66) is investigated 
at the LO of QCD. The quantitative role of NLO corrections and the implications of the (in)stability of 
existing calculations of (jfj^'-'i'tj bb) for explaining the observed puzzle is discussed in Section 5. The 
conclusions are drawn in Section 6. 



2 Basic facts and formulae 

The basic quantity of perturbative QCD calculations, the renormalized color coupling as{fJ-), depends on the 
renormalization scale /i in a way governed by the equation 

where for n/ massless quarks (3o = 11 — 2n//3 and /3i = 102 — 38n//3. Its solutions depend beside also on 
the renormalization scheme (RS). At the NLO this RS can be specified via the parameter Ars corresponding 
to the renormalization scale for which as diverges. The coupling as(M) then solves the equation 



At the NLO the coupling is a function of the ratio /i/A^g and the variation of the RS for fixed scale /i is 
therefore equivalent to the variation of /i for fixed RS. To vary both the renormalization scale and scheme 
is legitimate, but redundant. Throughout the paper I will work in the conventional MS RS and vary the 
renormalization scale only. As we shall investigate the QCD predictions down to quite small values of the 
renormalization scale ^, the equation |5| will be solved numerically, rather than expanding its solution in 
inverse powers of ln(/x/A). 

The main difference between hard collisions of hadrons and photons comes from the fact that quark and 
gluon distribution functions of the photon 

"/ "/ 
I](x,M) = ^fe(x,M) + g,(a;,M)) , q^s{x,M) = ^ (e? - (e^)) {q,{x,M) +q^{x,M)) (3) 
1=1 1=1 
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Figure 2: Diagrams defining the pointlike parts of the nonsinglet quark and gluon distribution functions of 
the photon, represented by the symbols on the left. 



satisfy the system of coupled inhomogeneous evolution equations 



= 5^kq{M)+Pqq{M)®i:{M)+PqG{M)®G{M), 



<m{M) 
dh^ 

= kG{M) + PGq{M)®^{M)+PGG{M)®G{M), 

'^^^ = 5nsk,{M)+P^s{M)®q^s{M), 



where (5ns = 6n/ ((e"*) — (e^)^), (5s — 6nf{e'^) and 



kg{x,M) = — 



kG{x,M) = £ 



P^jix,M) = 



27r 



2n 



(4) 
(5) 
(6) 

(7) 
(8) 
(9) 



The lowest order inhomogeneous splitting function kq^\x) = {x'^ + (1 — x)'^) as well as the homogeneous 
splitting functions P^p (x) are unique, whereas all higher order splitting functions fcg-'' , /c^^ , pjf^ , j > 1 depend 
on the choice of the factorization scheme (FS). Although potentially important, I will not exploit this freedom 
and throughout this paper will stay within the MS FS. The equations (|4I6|I can be recast into evolution 



equations for qi{x,M),q^{x,M) and G{x,M) with inhomogeneous splitting functions k^q^' 



Due to the presence of the inhomogeneous terms on the r.h.s. of H4I()|) their general solutions can be written 
as a sum of a particular solution of the full inhomogeneous equations and a general solution, called hadron- 
like (HAD), of the corresponding homogeneous ones. A subset of the former resulting from the resummation 
of contributions of diagrams in Fig. |21 describing multiple parton emissions off the primary QED vertex 
7 ^ 99 and vanishing at M = AIq, are called point-like (PL) solutions. Due to the arbitrariness in the choice 
of Mq the separation 

D{x, M) = D^^{x, M, Mo) + D^^^{x, M, Mq). (10) 
is, however, ambiguous. The explicit form of the pointlike contribution to nonsinglet distribution function 



q^^{n,Mo,M) 



in 



as{M) 



1 - 



s{M) \ 



l-2pW(n)//3o- 



asiMo)J 



fc^g(n) 



2^/3ol-2P,(°^(n)//3o 



(11) 



is often claimed to imply that it behaves as 0{a/ag). However, the fact that as{M) appears in the denomi- 
nator of (7^g cannot be interpreted in this way 19' because switching QCD off by sending Ars for fixed 
M, Mq reduces, as expected, the expression (|ll|l to the purely QED contribution, corresponding to the first 
diagram in in the upper part of Fig. El 



q^^(x,M,Mo)-^ — k^°l{x)\n'-^ 



2n 



(12) 
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The form merely implies that for asymptotically large M the pointlike part g^g is proportional to InM. 

As emphasized long time ago by Politzer |2(JI there is no compelling reason for identifying the renormalization 
and factorization scales /x and M and one should therefore keep these scale as independent free parameters 
of any finite order perturbative calculations. 



3 QQ production in 77 collisions 

We shall first recall the general form of the perturbative expansion of the cross section 0-404(77 ~* QQ'^ ^) ^^'^ 
then discuss in detail the renormalization and factorization scale dependence of finite order approximations 
to the three QCD contributions to this cross section. 



3.1 General form of atot{ll QQ) 

In the calculations of refs. [ZllHlinij performed with fixed pole quark masses, the NLO QCD approximation 
to utotill QQ) is defined by taking into account the first two terms in the expansions of direct, as well 
as single and double resolved photon contributions. Up to the order and suppressing the dependence on 
77 collision energy W these expansions read 

^dir (M) = a^^^ + <ji2 a, (m) + 4?r (M, f^)al (m) + (M, ^^)al{fi) + ■ ■ ■ , (13) 
a,,(M) = aW(M)a,(M)+a(2)(A/,M)a2(/.)+a(3)(M,Ai)«3(^) + ..., (14) 

^dr(M) = <j^^{M)al{fi)+<jl^^{M,fi)aUf^) + ---. (15) 

Starting at order the direct photon contribution depends also on the factorization scale and therefore 
mixes with the single and double resolved photon ones. The first two terms in p3|l are, however, totally 
unrelated to any terms in (|14|) or H15|l . 

The approximations employed in 0015] include all terms that are currently known, so we cannot presently 
do better. On the other hand we should be aware of its theoretical deficiency. The fact that the first two 
terms of (|13I15|) start and end at different powers of as is usually justified by claiming that PDF of the 
photon, which appear in expressions for Csr' (M) and a^^' {M), behave as a /as- Consequently, the first 
terms in all three expressions (|13I15|I are claimed to be of order (as)" = 1 and the second ones of order 
as- However, as emphasized above and argued in detail in the term InM^ characterizing the large M 
behaviour of PDF of the photon comes from integration over the transverse degree of freedom of the purely 
QED vertex ^ ^ qq and cannot therefore be interpreted as 1/asiM). 



3.2 Direct photon contribution 

For proper treatment of the direct photon contribution (|13f) . the total cross section of e+e~ annihilations 
into hadrons at center-of-mass energy y/^ provides a suitable guidance. For nj massless quarks we have 

ahad(V^) - <yZ{^) + «.(M)4ad(^) + «'(M)4ad(V^/M) + • • ' = '^("^,(1 + r{VS)), (16) 

where the term a^^^ = (Aira'^/S) ^/ comes, similarly as a^l in H13|) . from pure QED, whereas genuine 

QCD effects are contained in the quantity 

r{^)^^f^^^\l + as{^i)n{Vs/^l) + ■■]. (17) 

TT L J 

For the purpose of QCD analysis of the quantity Ijltill it is a generally accepted practice to discard the lowest 
order term cr|^°|j and denote as the "leading order" the second term in H16|) . i.e. ChadO^s/Tr. The adjectives 
"LO" and "NLO" are thus reserved for genuine QCD effects described by r{\/S). The rationale for this 
terminology is simple: to work in a well-defined renormalization scheme of as requires including in (|17|l at 
least first two consecutive powers of as{fi)- The expHcit /x-dependence of ri{\/S / cancels to the order a1 
the implicit /x-dependence of the leading order term Ug (fJ.) /n in H17|) and thus guarantees that the derivative 
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Figure 3: Examples of diagrams describing the direct photon contribution to cr (77 QQ) up to the order 
a^Og. The soUd (dashed) Hnes denote hght (heavy) quarks. 

with respect to In /i of the sum of first two terms in (|17|l behaves as . For purely perturbative quantities like 
(|16() the association of the term "NLO QCD approximation" with a well-defined renormalization scheme is a 
generally accepted convention, worth retaining for any physical quantity, like the direct photon contribution 
CTdir in H13|l . Contrary to this practice, the calculation in refs. X«S«»2. consider the purely QED contribution 



1 



In 



1 + 

1-/3 



/3 1 + 



4m^ 



127re|c^ 
1^2 ' 



(18) 



where /3 — — Am^jW'^ , as the LO approximation. This is legitimate but implies that their NLO ap- 
proximation, includes only the lowest order term in ois and cannot therefore be associated to a well-defined 
renormalization scheme of as- For QCD analysis of (Tdir in a well-defined renormalization scheme the incor- 
poration of the third term in (|13|l . proportional to c?o?g^ is indispensable. 

At the order o?c?g the diagrams with light quarks appear and we can distinguish three classes of direct 
photon contributions, differing by the overall heavy quark charge factor CF: 

Class A: CF — eg. Comes from diagrams, like those in Fig. Ot-g, in which both photons couple to heavy 
QQ pairs. Despite the presence of mass singularities in contributions of individual diagrams coming 
from gluons and light quarks in the final state and from loops, the KLN theorem implies that at each 
order of as the sum of all contributions of this class to CTdir is finite. Note that the first as well as 
the second terms in Ijlbf) are also proportional to Cq and it is therefore this class of direct photon 
contributions that is needed for the calculation of CTdir to be performed in a well-defined RS. 

Class B: CF = eg. Comes from diagrams, like that in Fig. I^h, in which one of the photons couples to a 
heavy QQ and the other to a light qq pair. For massless light quarks this diagram has initial state mass 
singularity, which is removed by introducing the concept of the light quark (and gluon) distribution 
functions of the photon. The factorization scale dependence of the contribution of this diagram is then 
related to that of single resolved photon diagrams in Fig. 2Jl,c. 

Class C: CF — 1. Comes from diagrams in which both photons couple to light qq pairs, as those in Fig. 
In this case the analogous subtraction procedure relates it to the single resolved photon contribution 
of the diagram in Fig. 2f and double resolved photon contribution of the diagram in Fig. The 
classes B and C are thus needed to guarantee the factorization scale (and scheme) invariance of the 
single and double resolved photon contributions to order cP'c?!.. 

Because of different charge factors CF ^ the classes A, B and C do not mix under renormalization of as and 
factorization of mass singularities. As the diagrams in Fig. |3l3 and 01 give the same final state qqQQ-, we 
should consider their interference term as well, but it turns out that it does not contribute to the total cross 
section atotill QQ)- 



3.3 Resolved photon contribution 

The classes B and C of direct photon contributions of the order a^a^ are indispensable to render the sum of 
direct and resolved photon contributions factorization scale invariant up to order a^aj.. To see this in detail. 
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Figure 4: Examples of resolved photon diagrams involving the pointlike parts of PDF of the photon and the 
related direct photon diagrams. 



let us write the sum of first two terms in H14I15|I explicitly in terms of PDF and parton level cross sections 
2aliii) jj dxAyY,q,{x,M)q^{y,M) \af^ {xy) + a,{^l)a'^^ {xy,M, ^i) 

i 

al[^Ji) jj dxAyG{x,M)G{y,M) [a%[xy) + a,{^l)af(,{xy, M, ti) 

2al{y) j I dxAy^x, M)G{y, M)a^^^ {xy, M) (19) 

where qi runs over Uf quark flavors and the factors of two and four reflect the identity of beam particles 
and equality of contributions from quarks and antiquarks. Recalling the general form of the derivative 



dlnAf2 



AxWo{x,M) + dx 



dccdy 



J2 M)W^^ {x, M) + G(x, M)Wg{x, M) 



G{x, M)Giy, M)WGG{xy, M) + ^ q,{x, M%{y, M)W,-,ixy, My 



i:{x,M)G{y,M)W,G{xy,M) 

using (|4I6|I and denoting as = as (/i) , / = d//d In Af ^ we find 

aa'i 1 



(20) 



Wo{x,M) 

Wg,ix,M) 

Wg{x,M) 

WGGix,M) 

W,g{x,M) 
WgGix,M) 



^ I i^'o (-)-S (-) + efcf (-) E ^) I + • • • 

^ 1 47ra(2) (x) +jdz [P(°) (z)a;2 (xz) + 3e?aA:(°) (z)ag) (xz) 
^(^2na'^l{x) + y'dzF(°^(z)a«(xz)| + • • ■ 
^ {7rag^(x) + jdzP^%iz)a^^lixz)^ +■■■ 
^ \2nall\x) + 2 /dzP(«)(z)a(y (a;z) j + • • • 



^ \ 2vrag(:r) + I dzlp^'J {z)al^ (xz) + P^^l{z)4l{xz) 



(21) 

(22) 
(23) 
(24) 
(25) 
(26) 



Only the lowest order terms on the r.h.s. of 121I26|I are written out explicitly. All integrals in 122I26II go 

formally from to 1, but threshold behaviour of cross sections (Jij{xz) restricts the region to xz > Aitiq/W'^ . 
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The factorization scale invariance of (|19|l requires that its variation with respect to In is of higher order 
in as than the approximation itself. There is no dispute that direct photon contributions of classes B and C 
are needed to guarantee this property. The question is which terms on the r.h.s. of (|20|l must vanish if the 
approximation is defined by H19|l . 

In the conventional approach both q{M) and G{M) are claimed to be of order a/ as and the approximation 
H19(l thus of the order a^as, implying that only terms up to this order must vanish in H2()|l. This in turn 
means that the functions H22I23|I must vanish to order a1 and (|24I26|I to order respectively, which, indeed, 
they do ^. The fact that the expression on the r.h.s. of H21(l does not vanish is of no concern in this approach 
as it is manifestly of the order aal and thus supposedly of higher order than (|19l) itself. 

If, on the other hand, we take into account that quark and gluon distribution functions of the photon 
behave as q{M),G{M) cx a, we see that Wo is of the same order a'^a^ as the products qiWq-,GWc and 
other integrands on the r.h.s. of H20() and must therefore also vanish for theoretical consistency of the 
approximation H19|) . This, in turn, necessitates the inclusion of class B direct photon contributions of the 
order aa^, like those in Fig. 0l3,g, which provide the M-dependent terms the derivative of which cancels 
the first term in (|20|l involving the integral over Wq. Note that Wg^ in H22|l receives contributions from the 
derivatives of both single and double resolved photon diagrams, proportional to a-yo and aqq, respectively. 
This fact reflects the mixing of single and double resolved photon contributions, which starts at the order 
a^al and is due to the presence of the inhomogeneous splitting terms in the evolution equations H4l(j|) . For 
theoretical consistency of the sum of direct and resolved photon contribution up to the order aa^ only the 
lowest order double resolved photon contribution must be included. 



4 66 production at LEP2 

We now turn to the phenomenological analysis of bb production at LEP2, where the incoming leptons act as 
sources of transverse and longitudinal virtual photons, described by the fluxes 

where stands for photon virtuality. Although the kinematic region of the LEP data includes photon 
virtualities up to moderate Q^, the cross section of the inclusive process 

e~ ^ e~ bb + anything (29) 

is dominated by the production of the bb pair in the collision of two quasireal photons with very small Q^, 
typically (Q^) ~ 0.01 GeV^. For such small the cross sections of hard processes involving longitudinal 
virtual photons, which are proportional to Q^, are expected to be negligible compared to those of transverse 
virtual photons. When talking about the production of bb in e+e~ collisions we shall always mean in 
association with the e^e~ pair, but for brevity of notation shall drop this latter specification, writing 
o'tot(e^e~ bb) instead of crtot(e+e^ e+e~66). 

Although the data are available only for cross sections integrated over the whole phase space, we shall discuss 
the contributions dcrfc(e+e~ bb)/dW of individual processes as functions of 77 collision energy W. The 
shapes of these contributions can altenatively be characterized by the functions 

^ dak{e+e-^bb) r< - a dCTfc(e+c- ^ 66) 



Fk{W)= dw ^ GkiW)= dw ^ (30) 



2m, Jw dw 



which quantify how much of a given contribution is located in the region up to a given W {Fk(W)) or 
above it {Gk{W)). As the available data are not copious enough to measure the differential distribution 
dcr(e+e~ 66)/dW^ the theoretical analysis of the distributions (|30|l might allow us to invent a strategy how 

^The latter condition is actually the same as for QQ production in hadron-liadron collisions. 
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Figure 5: Left: The distributions dak/dW corresponding to the pure QED contribution together with three 
lowest order QCD contributions: single resolved (SR), double resolved (DR), and direct (DIR). The sum of 
all four contributions is shown as the upper dotted curve. Right: the comparison of pure QED contribution 
with the sum of three lowest order QCD ones. All curves were obtained for nib — 4.75 GeV, ^/S = 200 GeV 
and Q2^„^ ^ 4 GeV^ using GRV LO PDF of the photon and setting A^") = 0.27 GeV. 



to separate the kinematic region of accessible W into two parts, each dominated by a particular contribution. 
The relative importance of the individual contributions as a function of W is determined by the ratia 



TkiW) 



dcrfc(e+e- 



dW 



dgtof(e+e bb) 
dW 



(31) 



4.1 QED contribution 



The pure QED contribution to crfot(e"''e 



IS given as 



do-QED(' 



dW 



Qa^el A{W) 



where 

A{W) 



Ami 



8m* 



dydz6 ( — - 



1 + (1 - y)2 



In 



l + (3 



/3 1 + 



m'^y^ 



In- 



)2 

C max 



(l-z) 



(32) 



(33) 



results from convolution of of photon fluxes 127|) , integrated over the virtualities up to Q^^ix ■ The convolution 
H33|) can easily be performed analytically and the result inserted into H32I) . In Fig. [Slwe display by the solid 
curve the result of evaluating for mb = 4.75 GeV, = 200 GeV and Ql^^x = 4 GeV. The distribution 
vanishes at the threshold W = 2mb due to the threshold behaviour of the cross section ((TH)l . peaks at about 
W = 12 GeV and then drops rapidly off due to the fast decrease of both the photon flux (|27l) and (|18|l . 
Integrating the distributions in Fig. yields the values in the fourth column of Tabled 



4.2 Leading order QCD corrections 

QCD corrections to pure QED expression (|18|l are of three types: direct (dir), single resolved (sr) and double 
resolved (dr). We shall first discuss the lowest order contributions to all three types of QCD corrections. As 
in the case of pure QED contribution, these corrections are given as convolutions of the photon flux H27() with 
the appropriate partonic cross sections. In all calculations u, d, s and c quarks were considered as intrinsic 
in the photon and nf = A was taken in the expression for as{fJ-)- 
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Parameters 


QED 


LO QCD 


Total 


rub AW PDF 




DIR 


SR 


DR 




4.75 0.27 GRV LO 


1.27 


0.473 


1.415 


0.121 


3.28 


4.5 0.27 GRV LO 


1.40 


0.478 


1.746 


0.146 


3.77 


4.75 0.35 GRV LO 


1.27 


0.520 


1.542 


0.141 


3.47 


4.75 0.27 SASID 


1.27 


0.473 


0.904 


0.077 


2.73 



Table 1: The integrated cross sections cr(e+e~ bb,S) for ^/S = 200 GeV and Q^ax — 4 GeV^, corre- 
sponding to the distributions in Fig. El The renormahzation and factorization scales /i and M we identified 
and set equal to rrib. LO form of was used. All cross sections are in picobarns. 



4.2.1 Direct photon contribution 

The W dependence of the leading order QCD correction is given as the product 

= W™.) (34) 

of the convolution ^(W^) of photon fluxes and the lowest order QCD contribution as{n)(J^^^^{W). At this 
order the direct photon contribution cr^|^Q:s(/i), which comes form real or virtual emission of one gluon, is 

exclusively of class A. The function crjj^(VK/mf,) has been calculated in, for instance, As it is just the 
first term in the series in positive powers of as (/x) , the value of the renormahzation scale jj, in the argument 
of as{p) is completely arbitrary. The resulting VF-dependence, evaluated for /i — mi, and shown in Fig. [SJ 
is peaked even more sharply at small W than the pure QED contribution (|32|l . This reflects the fact that 
the cross section o"^j^^(VK/mb) does not vanish at the threshold W ~ 2mf, as does a^^^^{Wlmi,). 



4.2.2 Resolved photon contribution 

The leading-order single and double resolved photon contributions, were computed with HERWIG Monte 
Carlo event generator, which implements the appropriate LO cross sections of the processes 

7 + G ^ h + b, (35) 
G + G b + b, q + q^b + b, (36) 

where q = u,d, s,c stand for intrinsic quarks in the photon, and convolutes them with photon fluxes and PDF 
of the quasireal photon(s). In HERWIG the renormahzation and factorization scales fi and AI are identified 
and set equal to an expression which is approximately equal the transverse mass jj. = M — Mt = \/ + AP. 
In LEP2 energy range the mean (Mj) depends weakly on W with, approximately, (Mj) ~ 7 GeV. 

Results of the calculations in which the LO GRV PDF of the photon, the LO expression for as{^) with 
A(4) = 0.27 GeV and mb = 4.75 GeV were used, are shown in Fig. As expected, the corresponding 
distributions are much broader than those of pure QED or LO QCD direct contributions. 



4.3 Comparison of individual contributions 

The comparison of the distributions da/c/dM^, i^fe(VF) and Gk{W), corresponding to four individual contri- 
butions, displayed in Figs. (SjandEland summarized in Tabled reveals large difference in their shapes and 
magnitude. Specifically we conclude that 

• The pure QED as well as the LO direct photon contributions peak at very small W and are basically 
negligible above ~ 50 GeV. For instance, the left plot of Fig. shows that 95% of the QED 
contribution comes from the region < 30 GeV. 

• The onset of single as well as double resolved photon contributions is much slower, but these distribu- 
tions are, on the other hand, markedly broader. 
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Figure 6: Left: solid (dashed) curves show the partially integrated cross sections Fk{W) {Gk{W)) defined 
in l|30() for QED and three LO QCD contributions. Right: The relative contributions rk{W) defined in H31() 
for the same four contributions. 



• The double resolved photon contribution is negligible everywhere. 

• The pure QED and single resolved photon contributions are of comparable size and together make up 
about 85% of the total integrated cross section, 

• up to about ~ 30 GeV, datot/<iW is dominated by pure QED contribution, whereas for W > 30 
GeV, QCD contributions take over. 

The numbers given in Table 1 correspond to standard colored quarks with fractional electric charges. In 
the excess of data over standard theoretical calculations is interpreted as evidence for Hahn-Nambu integer 
quark charges. Applied to the case of b-quark, the author of argues that the correct way of calculating 
the charge factor in H32|l is not the usual 3e^ = 1/27, but ^b'')*/'^ = 1/3 = 9/27, where the sum runs 

over the three Hahn-Nambu integer b-quark charges e^'\ which are 0, 0, —1 respectively. The results is thus 
9 times bigger than that of the standard calculation. I think his argument for first summing over the quark 
colours and then taking the fourth power is wrong ^, but I mention it here as illustration of the merit of 
separating the data into at least two regions of W. Were the author of jJHI right, the whole discrepancy 
would have come from the region of small W, where QED contribution dominates. 

On the other hand, were the light gluino production responsible for the observed excess, the latter would 
have to come from the region of W dominated by the double resolved photon contribution. Although the 
energy dependence of the gluon-gluon fusion to gluino- antigluino may be slightly different than those of 
G + G ^ QQ or qq — > QQ, it is clear that the basic shape of the ly-distribution is given by the convolution 
of the photon fluxes H27I28|I and the gluon distribution function of the photon, which are the same in both 
types of processes. 

The above observations underline the fact that in order to pin down the possible origins of the excess of the 
integrated cross section crtot(e+e^ e'^e~bb) over theoretical calculations, it would be very helpful if the 
data could be separated at least into two subsamples according to their hadronic energy W, say W < 30 
GeV and > 30 GeV. 

The magnitude of the contributions discussed in the preceding subsection depend, beside the e+e~ cms 
energy v^, on a number of input parameters: the numerical values of m?,, Aqq^, Q^axj the selection of PDF 
and the choice of the renormalization and factorization scales fi and M. In all the calculation reported above 
we set ji = M = rrit- The central calculation was performed for y/S — 200 GeV, Q^ax — 4 GeV^, mt = 4.75 

^The correct value of the charge factor in t'i'2l for the Hahn-Nambu integer charge b-quark equals 1 and would thus yield 
even bigger enhancement than that suggested in llbl . 
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GeV, A('*) = 0.27 GeV using the GRV LO PDF of the photon. To see the sensitivity of the LO resuhs to 
these assumptions we varied some of these parameters: 

• rrifc was lowered to irib = 4.5 GeV, 
. A(4) was increased to 0.35 GeV, 

• GRV set of PDF of the photon was replaced with that of Schuler-Sjostrand set SASID. 

The choice of Q^ax — 4 GeV^ corresponds roughly to the usual cuts imposed on the LEP2 data and could 
therefore be also adjusted to specific conditions of a given experiment. 

The results of the calculations of tTtot(e+e~ e^e'bb), corresponding to different sets of input parameters 
specified above, are listed in Table Lowering nib increases all four contributions, as does, except for the 
pure QED one, increasing A'"'). SASID PDF yield markedly lower results for single and double resolved 
photon contributions. It is, however, clear that varying the input parameters within reasonable bounds does 
not bring the sum of lowest order QED and QCD calculations significantly closer to the data. 



5 Can the NLO QCD corrections solve the puzzle? 

With the sum of lowest order QED and QCD contributions to a^§{e'^e~ e+e^65) way below the data we 
shall now address the question whether the next-to-leading order QCD corrections can at least partly bridge 
the gap between data a theory. 



5.1 Direct photon contribution 

The the sum of the second and third terms in (|13|) 
the ratio W/rub, as 



can be written, suppressing the dependence of cr^*j|, on 



^NLO 
Mir 



(1) 



i^i) [l + rl{n/mb)as{^i)] 



(37) 



Note that W'^o'^lr ^ unique function of the ratio W/rub and the NLO coefhcient ri{W, nib, m) can be written 
as a function of W/rub and nib/fJ-. The first term in 137|l is a monotonous function of the renormalization 
scale /i, spanning the whole interval between zero and infinity. As emphasized in Section 13.21 one needs to 
include at least the term a^cr^'^j'' to make the expression (IS7|l of genuine next-to- leading order. The class A 
of order a^a^ direct photon contributions is needed for this purpose. 

The renormalization scale invariance of cr^^^ implies the following general form of ri : 



ri(W^/mf,,^/mf„RS) = 7^ In -rj p{W/mb), 

Aj^s 



(38) 



where p{W/mb) is a renormalization scale and scheme invariant |17| . which, however, depends beside the 
ratio W/irib also on the numerical value of the ratio m;,/A'^'*^. It can be evaluated using the results of a 
calculation in any renormalization scheme RS 



p{W/mb, mb/A^) ^ ^ In i^l^ - n [W/mb, 1, RS) 

^RS 



(39) 



and its numerical value governs basic features of the scale dependence of H37|l : 



p > 0: the NLO approximation (|37|l considered as a function of p exhibits a local maximum, where 
dcr^;^°/d/j, = and where the prediction is thus most stable. This point, preferred by the Principle of 
Minimal Sensitivity J7], is also very close to the point for which ri = 0, which is selected by method of 
Effective Charges The value of <J^^'~' at this point is proportional to l/p implying very large NLO 
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Figure 7: a) The renormalization scale dependence of the leading (solid curve) and next-to-leading order 
contributions to the generic quantity (|42|l for different values of p. The dashed curves correspond to p > 0, 
the dotted ones to p < 0, the dash-dotted to p = 0; b) the same as in a) but plotted in a linear scale of /z; 
c) graphical representation of H42|) in three different renormalazation schemes and for p > 0; d) the shape of 
the NLO expression (|42() as a function of as- 



corrections for small p. Inserting the appropriate numbers for Uf = A, mi, = 4.75 GeV and A-|^^ = 0.27 
GeV, we get 

p{W/mb) ^3.88- ri(W/mb,l,MS). (40) 
The coefficient ri(W/mb, 1,MS) thus does not have to be large to get small, or even negative pi 

p < 0: CTjj^'-* becomes a monotonous function of /i, similarly to cr^P. In fact, taking the derivative of cr^^'-' 
with respect to p one finds that for p < it is actually even steeper than a^P , given by the first term in 
(IS7|) . Consequently, for negative p going to the NLO does not improve the stability of the calculation, 
but quite on the contrary! 

The above features are straightforward to see assuming /3i 0. This assumption simplifies the relevant 
formulae, but nothing essential depends on it. Setting c = in |(2Jl allows us to write explicitly 

An 

- /?oln(pVAL) ^''^ 
which, inserting this expression into (|37|l and taking into account (jJSJ), gives 



„NLO 
C^dir 



= a^;,^a.(p)[2-pa,(p)], (42) 



In Fig. [3 we plot the dependence of the generic NLO quantity (|42|l for several values, positive as well as 
negative, of p and in different renormalization schemes. Several conclusions can be drawn from this figure: ^ 

• For p < cr^ir ^ steeper function of p than a^P . For positive p, on the other hand, (|^ exhibits a 
local maximum at a™°^ = 1/p. 



•^AU conclusions mentioned below are well known in the context of perturbative quantities of the form 1871 depending on 
the renormalization scale only. We recall them because they will be used in the analysis of the resolved photon contribution, 
where the interplay between the renormalization and factorization scales complicates the situation. 
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• For negative as well as positive values of p, t^j]:''^ cx 1/ In as — s- oo. For negative p this implies that 
there is no region of local stability. However, when plotted on a linear scale of /i in a limited interval 
the weak logarithmic does, as illustrated in Fig. [Tjp, fake the local quasistability. 

• The curve representing cr^^'^ depends, as shown in Fig. [7|: for p > 0, on the chosen renormalization 
scheme. Setting p equal to some "natural" physical scale Q therefore does not resolve the renormal- 
ization scale ambiguity as in different RS we get different results. However, although also the position 
of the local maximum depends on the choice of the RS, the value of the NLO approximation H42(l at 
this maximum does not! The same holds for the intersection of the LO and NLO curves |l21j . 

• Instead of varying both the renormalization scale and scheme, which is legitimate but redundant, we 
may use the couplant as itself for labeling the different predictions of (|42|l . Instead of an infinite set 
of curves describing the ^-dependence of o'j;^'-*(//, RS) in different RS, we get for each p a single curve 
displayed in Fig. |7Ji. 

We wish to emphasize that whereas there are natural physical scales in hard collisions, there is nothing like 
the "natural" renormalization scheme. The fact that the choice of renormalization scheme is in principle 
as important as that of the renormalization scale is often not fully appreciated. For instance, the standard 
way of estimating the importance of higher order corrections beyond the NLO approximation employs the 
ratio (called "k- factor") k =NLO/LO evaluated for some "natural" renormalization scale. If this ratio is 
significantly larger than 1, perturbation theory is deemed unreliable. However, the "k- factors" themselves 
depend, for a chosen renormalization scale, on the choice of the renormalization scheme and so does therefore 
also the importance of higher order corrections! The usual procedure of estimating the convergence of a 
given perturbation expansion by the value of the corresponding "k- factor" evaluated for some "natural" 
renormalization scale in the calculationally convenient MS renormalization scheme is thus entirely ad hoc. 

As the term a^^^. in H37|l has not yet been calculated, we cannot associate the class A direct photon contribu- 
tion, given by the first term in H37() . to a well-defined renormalization scheme. Moreover, as the magnitude 
of the NLO corrections in H37|) is determined by the ratio o-l^^iW/mh,l)/al^^{W/'mh) of two functions of 
W/nib, which may depend on W/irib in different ways, the coefficient ri may be very large even when both 
the numerator and denominator are on average of comparable and small magnitude. The size of NLO correc- 
tion may also depend sensitively on the ratio W/nib- All this indicates that without the knowledge of class 
A direct photon contribution of the order a^a^, we cannot make a meaningful estimate of the importance 
of higher order corrections to H37() . 



5.2 Resolved photon contribution 

As in LEP2 energy range the double resolved photon one is numerically negligible, only the single resolved 
photon contribution will discussed in detail below. As shown in Fig. [SJ the spectrum of the contributions 
dasr/dW peaks at about W = 30 GeV with the mean value {W) = 65 GeV. The properties of the measured 
cross section (T^'"'-'(e+e^ — + e+e~66) will therefore be determined primarily by those of cr^'"° (77 — > bb) in 
the energy range 30 < < 65 GeV. All the results shown below for the single resolved photon contribution 

aNLO(H/^ ^) ^ 2as{p) J dxG{x, M) [cj%{x) + as{p)a%(x, M, p) 

Aalip) I dx^(Z,(a;,M)42)(x,M) (43) 

are based on the formulae for the partonic cross sections ajj^-* as given in Even if the reader does 

not agree with our claim that the approximations employed in |7||S1|5] do not constitute complete NLO 
approximation, it is certainly important to understand quantitatively their renormalization and factorization 
scale dependence. Because the expressions for ct^q as given in correspond to p = Af, we have restored 

its separate dependence on p and M by adding to (t^^q{x, M, M) the term (/3o/47r) cr^g, ln(/Li^/M^). Note that 
for each value of M the expression 143(1 has, as far as the /i-dependence is concerned, the form of the NLO 
expression (pTfjl . In all calculations the GRV HO set of PDF of the photon and A^ = 0.27 GeV were used. 
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Figure 8: a) The scale dependence of the conventional NLO approximation af^^'~' {W, M , M) for W = 40 
GeV (solid curve) together with the contributions of the 7G (dashed) and ^q/q (dotted) channels. The LO 
approximation is shown for comparison by the dash-dotted curve, b) a^^^'^ {W, M , kM) (solid curve) and 
al'PiW, M, kM) (dash-dotted) for W = AO GeV and three values of k = 0.5, 1.2. 
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Figure 9: The sur face and contour plots of <^"(VF, M, fi) for W = 40 GeV 



We first follow the conventional procedure and set M = ^. The resulting (common) scale dependence of 
the expression (|43|l . together with those of the quark and gluon contributions to it, are shown in Fig. |St 
for W ^ 40 GcV. Overlaid for comparison is also the LO approximation, given by the first term in (|43|l . 
We note the different scale dependence of the and 75 channels, the latter turning negative for Af > 6 
GeV, but the most important observation concerns the fact that the conventional NLO approximation H43|l 
is a monotonously decreasing function of the scale. Moreover, it falls off even more steeply than the LO 
expression! In other words in going from the leading to the next-to-leading order the sensitivity to the 
scale variation increases, rather than decreases, as one might expect (and hope)! Recalling the discussion 
in Section 5.1 one should, however, not be surprised. To check how much this feature depends on setting 
exactly /i = M, we plot in Fig. the scale dependence of (t^.'"^(W, M, /i = kM) for standard choices of 
K = 0.5, 1, 2. Clearly, the above conclusion is independent of k in this range. 

The steep and monotonous scale dependence of (T^^°(W, M, 11 = kM) is a clear warning that the conventional 
NLO approximation is highly unstable. To see what happens if we relax the usual but arbitrary identification 
^ = kM we plot in Fig. |51 the surface and contour plots representing the full M and fi dependence of 
CT^.^° ( W, M, /x) as given in eq. Contrary to analogous process in antiproton-proton coUisions it 

does not exhibit a saddle point, where the derivatives with respect to both M and /i would vanish, but Fig. 
El seems to indicate some sort of stability region at large scales, say for M > 10 GeV, /i > 20 GeV. This 



14 




Figure 10: The renormalization scale dependence of afj^'~'{W, M, /x) for fixed values of the factorization scale 
(a) and vice versa: the factorization scale dependence of af^^*^ {W, M , ii) for the same set of fixed values of 
the renormalization scale /i (b). All calculations correspond to = 40 GeV. In a) the ordering from above 
of the curves at ^ = 30 GeV corresponds to M = 30, 16, 10, 7, 4.75, 3, 2, 1 GeV, in b) the curves correspond 
at M = 0.7 GeV to the same sequence from below. 

impression is, however, misleading as becomes clear if we plot in Fig. ^|the slices of the surface plot in Fig. 
IHK along both axis and recall the discussion of Section 5.1. For each fixed value of the factorization scale 
M the expression (I43|l has a form of the NLO expression as far as the renormalization scale /i is concerned. 
Comparing the curves in Fig. llOh with those ofE^L we see that for M < 4.2 GeV af^^'~' {M , fi) corresponds 
to negative p in H42|l and thus exhibits no local stability point. For higher M the local maximum in the 
^-dependence of (T^^'-'(Af, fi) exists at the associated fi^axiM). The Af-dependence of cr^^'-'(M, Hmax{M)), 
shown in Fig. llOb . by the dotted curve, is, however, even steeper that those at fixed M. The above plots and 
conclusions concerned the results at one typical value of W, but their essence holds for the whole interval 
relevant for LEP2 data. 

We thus conclude that in the energy range relevant for LEP2 data the renormalization and factorization 
scale dependence of the conventional NLO calculations of single resolved photon contribution to the total 
cross section atot (77 bb) exhibits no stability region, either as a function of the common scale /i = kM or 
as fully two dimensional function of /i and M. 

6 Summary and Conclusions 

We have argued that in order to understand the origin of the discrepancy between LEP2 data on bb production 
in 77 collisions and the current theoretical calculations, two ingredients are needed. On the experimental 
side, the separation of data into at least two bins of the hadronic energy W, say W < 30 GeV and W >30 
GeV, could be instrumental in pinning down the possible mechanisms or phenomena responsible for the 
observed excess. 

On the theoretical side, the evaluation of the direct photon contribution of the order a^a^ is needed to make 
the existing theoretical expressions of genuine next-to-leading order in as- In their absence, the existing 
NLO QCD calculations are highly sensitive to the variation of renormalization and factorization scale and 
thus inherently unreliable. 
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Abstract: Understanding the data on the total cross section (Ttoi(e^e~ — >e"'"e~66) measured 
at LEP2 represents a serious challenge for perturbative QCD. In order to unravel the origins of 
the discrepancy between data and theory, we investigate the dependence of four contributions 
to this cross section on 77 collision energy. As the reliability of the existing calculations of 
atot{^^^~ — >e^e~66) depends, among other things, on the stability of calculations of the cross 
section atoti^^ bb) with respect to variations of the renormalization and factorization scales, 
we investigate this aspect in detail. We show that in most of the region relevant for the LEP2 
data the existing QCD calculations of C7tot(77 bb) do not exhibit a region of local stability 
and should thus be taken with caution. The source of this instability is suggested and its 
phenomenological implications for LEP2 data are discussed. 
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1. Introduction 

Heavy quark production in hard collisions of hadrons, leptons and photons has been considered 
as a clean test of perturbative QCD. It has therefore come as a surprise that the first data on 
the bb production in pp collisions at the Tevatron [1,2], 7p collisions at HERA [3,4] and 77 

collisions at LEP2 [5,6] have turned out to lie significantly and systematically above theoretical 
calculations. The disagreement between data [5, 6] and theory [7-9] was particularly puzzling 
for the collisions of two quasireal photons at LEP2. 

The arrival of new data on bb production in ep collisions at HERA [10], shown in the 
left part of Fig. 1 as solid squares, have further complicated the situation. In the range of 
moderate ~ 80 GeV^ the new ZEUS data [10] are in reasonable agreement with NLO QCD 
predictions and also in the photoproduction region the excess of the new data over theory is 
substantially smaller then that in the older data. As a result, there is now an inconsistency 
between new ZEUS and older HI results [3] for moderate Q^, but the situation remains unclear 
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Figure 1: The current situation with data on 56 production in ep and e+e cohisions, including the 
most recent data of ZEUS [10] and DELPHI [13]. 



also in the photoproduction region. For pp collisions the progress on the theoretical side [11,12] 
has significantly reduced the discrepancy observed at the Tevatron. 

On the other hand, the problem of understanding the bb production in 77 collisions remains. 
The preliminary DELPHI data presented this spring at PHOTON 2003 conference [13] and 
reproduced in Fig. 1, are in striking agreement with the older L3 and OPAL data. The central 
values of all three experiments are almost identical which strongly supports the reliability of these 
measurements. Contrary to the case of analogous discrepancy in antiproton-proton collisions at 
the Tevatron, there have been few suggestions how to explain the sizable excess of data over 
current theory in 77 collisions. Neither the use of unintegrated parton distribution functions [14], 
nor the production of supersymmetric particles [15], proposed for explaining an analogous excess 
in antiproton-proton collisions, are of much help for LEP2 data, primarily because of low 77 
energies involved. Quite recently, however, this discrepancy has been interpreted as an evidence 
for integer quark charges [16]. We will come back to this suggestion in Section 4.1. 

In [12] we have investigated the sensitivity of QCD calculations of atotiPP bb; M, //) to the 
variation of the renormalization and factorization scales /i and M. In particular we have argued 
that in order to arrive at locally stable results [17] these two scales must be kept independent. 
We have furthermore shown that in the Tevatron energy range the position of the saddle point 
of the cross section atot {pp bb; S, M, /i) lies far away from the "diagonal" = M used in all 
existing calculations. Using the NLO prediction at the saddle point instead of the conventional 
choice ^ = M = mb enhances the theoretical prediction in the Tevatron energy range by a factor 
of about 2, which may help explaining the excess of data over NLO QCD predictions. 

In this paper similar analysis is performed for 77 collisions at the total centre of mass 
energies W relevant for existing LEP2 data. The specific features of the theoretical description 
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of QQ production in 77 collisions have been discussed in [18]. However, as all three experiments 
at LEP2 have measured merely an integral over the cross section (7(77 bb, W) weighted by 
the product of photon fluxes inside the beam electrons and positrons, is is also important to 
understand the VF-dependence of the four individual contributions to it. 

The paper is organized as follows. The basic facts and formulae relevant for the quanti- 
tative investigation of rcnormalization and factorization scale dependence of finite order QCD 
approximations arc collected in Section 2. This is followed in Section 3 by the discussion of 
the general form of atotill bb; W, M , . In Section 4 the VF-dependence of the four con- 
tributions to the cross section crtotle+e" — >-e+e~66) is investigated at the LO of QCD. The 
quantitative role of NLO corrections and the implications of the (in)stability of existing cal- 
culations of Ctoi'*^(77 bb) for explaining the observed puzzle is discussed in Section 5. The 
conclusions are drawn in Section 6. 



2. Basic facts and formulae 

The basic quantity of pcrturbative QCD calculations, the renormalized color coupling as{n), 
depends on the rcnormalization scale in a way governed by the equation 

^ - = -^«^(.) - ^«3(^) + . . , (2.1) 

where for Uf massless quarks /3o = 11 — 2n//3 and Pi = 102 — 38ra//3. Its solutions depend 
beside /x also on the renormalization scheme (RS). At the NLO this RS can be specified via the 
parameter Ars corresponding to the renormalization scale for which ag diverges. The coupling 
as(/x) then solves the equation 

At the NLO the coupling is a function of the ratio ^/Ans and the variation of the RS 

for fixed scale ^ is therefore equivalent to the variation of /j, for fixed RS. To vary both the 

renormalization scale and scheme is legitimate, but redundant. Throughout the paper I will work 

in the conventional MS RS and vary the renormalization scale // only. As we shall investigate 

the QCD predictions down to quite small values of the renormalization scale ^, the equation 

(2.2) will be solved numerically, rather than expanding its solution in inverse powers of ln(/i/A). 

The main difference between hard collisions of hadrons and photons comes from the fact 

that quark and gluon distribution functions of the photon 

"/ "/ 

S(x, M) ^ Yl ^) + ^)) ' l^six, M) ^ - ^^')) (*(^' ^) + ^)) (2.3) 

i=l i=l 

satisfy the system of coupled inhomogeneous evolution equations 

= S^kqiM) + Pgg{M) S(M) + PqciM) G{M), (2.4) 

= kaiM) + PcgiM) E(M) + Pgg{M) G(M), (2.5) 

= ^mk,{M) + Pns(M) gNs(M), (2.6) 
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where (5ns = 6n/ ((e^) — (e^)^), 6^, = Qnf{e^) and 



kg{x,M) 
kG{x,M) 
Pij{x,M) 



a 

a 
2^ 



«s(M)ji) 



27r 



«.(M) 
27r 



(2.7) 
(2.8) 
(2.9) 



The lowest order inhomogcncous splitting function k, 



(0), 



= (x^ + (1 — x)^) as well as the 
homogeneous splitting functions P-j' (x) are unique, whereas all higher order splitting functions 

kq ^ , kQ^ , Pjj^^, J > 1 depend on the choice of the factorization scheme (FS) . Although potentially 
important, I will not exploit this freedom and throughout this paper will stay within the MS 
FS. The equations (2.4-2.6) can be recast into evolution equations for qi{x, M),q^{x, M) and 



G{x, M) with inhomogeneous splitting functions k^^ = 3efk, 

Due to the presence of the inhomogeneous terms on the r.h.s. of (2.4-2.6) their general 
solutions can be written as a sum of a particular solution of the full inhomogeneous equations 
and a general solution, called hadron-like (HAD), of the corresponding homogeneous ones. A 
subset of the former resulting from the rcsummation of contributions of diagrams in Fig. 2 
describing multiple parton emissions off the primary QED vertex 7 — > and vanishing at 
M = Mo, are called point-like (PL) solutions. Due to the arbitrariness in the choice of Mq the 



.21.(0) 




Figure 2: Diagrams defining the pointlike parts of the nonsinglet quark and gluon distribution functions 
of the photon, represented by the symbols on the left. 



separation 



D{x, M) = D^^{x, M, Mo) + D^^"{x, M, Mq). 



(2.10) 



is, however, ambiguous. The explicit form of the pointlike contribution to nonsinglet distribution 
function 

l-2P^°\n)/po' 



Ait 



as{M) 



1 



as{M) \ 



(2.11) 
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where 

«Ns(n) ^ (2.12) 

is often claimed to show that it behaves as 0{a/as)- However, the fact that as{M) appears in 
the denominator of g^g cannot be interpreted in this way [19] because switching QCD off by 
sending Ars for fixed M, Mq reduces, as expected, the expression (2.11) to the purely QED 
contribution, corresponding to the first diagram in in the upper part of Fig. 2 

ql\{x,M,Mo) ^^4°)(x)ln^. (2.13) 

The form (2.11) merely implies that for asymptotically large M the pointlike part g^g behaves 
as In M. 

As emphasized long time ago by Politzer [20] there is no compelling reason for identifying 
the renormalization and factorization scales jjl and M and one should therefore keep these scale 
as independent free parameters of any finite order perturbative calculations. 



3. QQ production in 77 collisions 

We shall first recall the general form of QCD expression for atoti.ll ~^ QQ) ^'^id then discuss in 
detail the renormalization and factorization scale dependence of its finite order approximations. 

3.1 General form of atotill QQ) 

In the calculations of refs. [7-9], performed with fixed pole quark masses, the NLO QCD approx- 
imation to (Ttotill QQ) is defined by taking into account the first two terms in the expansions 
of direct, as well as single and double resolved photon contributions 

^dir(M) = + a^^asifi) + a^^^ (M, Ai)a^(/x) + ajg (M, /z)«3(^) + . . . , (3.I) 
asr{M) = a^\M)as{fi) + a^\M,p,)a^,if,) + a^\M,fi)al{fi) + ---, (3.2) 
^dr(M) = c7i?(M)a2(^)+^(3)(M,/x)af(/x) + ---. (3.3) 

Starting at order the direct photon contribution depends also on the factorization scale and 
therefore mixes with the single and double resolved photon ones. The first two terms in (3.1) 
are, however, totally unrelated to any terms in (3.2) or (3.3). 

The approximations employed in [7-9] include all terms that are currently known, so we 
cannot do better at this moment. On the other hand wc should be aware of its theoretical 
deficiency. The fact that the first two terms of (3.1-3.3) start and end at different powers of ag 

(1 2) 

is usually justified by claiming that PDF of the photon, which appear in expressions for Usr' {M) 
and a'^'^\M), behave as a/ag. Consequently, the first terms in all three expressions (3.1-3.3) 
are claimed to be of order (a^)*^ = 1 and the second ones of order as- However, as emphasized 
above and argued in detail in [19], the term InM^ characterizing the large M behaviour of PDF 
of the photon comes from integration over the transverse degree of freedom of the purely QED 
vertex 7 — > and cannot therefore be interpreted as l/a^ (M) . 
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3.2 Direct photon contribution 

For proper treatment of the direct photon contribution (3.1) to (Jtot{ll QQ)-, the total cross 
section of e+e" annihilations into hadrons at center-of-mass energy \fS provides a suitable 
guidance. For Uf massless quarks we have 

CThad(V^) = 4ad(V^) + + «^()«)4?d(^M) + • • • = 4ad(l + (3-4) 



(47ra^/3S') Cj comes, similarly as u^j] in (3.1), from pure 



where the term fij^^^ 

(^) 

QED, whereas genuine QCD effects are contained in the quantity r(\/5) 



asip.) 



TT 



l + as{fJ,)rl{^/S/f^) + ■■■ 



(3.5) 



For the purpose of QCD analysis of the quantity (3.4) it is a generally accepted practice to 
discard the lowest order term cr^^^ and denote as the "leading order" the second term in (3.4), 
i.e. <7^ad'^*(^)/^' "^^^ adjectives "LO" and "NLO" are thus reserved for genuine QCD effects 
described by r{VS). The rationale for this terminology is simple: to work in a well-defined 
renormalization scheme of requires including in (3.5) at least first two consecutive powers of 
asiiJ,). The explicit )U-dependence of ri(\/S'//i) cancels to the order the implicit /x-dependence 
of the leading order term as(/u) /vr in (3.5) and thus guarantees that the derivative with respect to 
In ij, of the sum of first two terms in (3.5) behaves as a^. For purely pcrturbative quantities like 
(3.4) the association of the term "NLO QCD approximation" with a well-defined renormalization 
scheme is a generally accepted convention, worth retaining for any physical quantity, like the 
direct photon contribution adir in (3.1). 



L, L, L, U ^ 



r^-- ;^ T 

J"""""" r^"""'*' r^"""*"" i^nnnnr wr\^ 5! 

j------ j--""-- '^"^^ 




d) e) f) g 



Figure 3: Examples of diagrams describing the direct photon contribution to (7(77 QQ) up to the 
order a^a^. The solid (dashed) lines denote light (heavy) quarks. 



Contrary to this practice, the calculation in refs. [7-9] consider the purely QED contribution 



1 + 



4mi 



8mt 



In 



1 + /? 
1-/3 



/3 1 + 



4m^ 



(To = 



(3.6) 



where /? = y 1 — 4m^/H^^, as the LO approximation. This is legitimate but implies that their 
NLO approximation, includes only the lowest order term in and cannot therefore be associated 
to a well-defined renormalization scheme of ag even if the NLO expression for ag is used therein. 
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For QCD analysis of ddir in a well-defined renormalization scheme the incorporation of the third 
term in (3.1), proportional to a^a^, is indispensable. 

At the order a^a^ the diagrams with light quarks appear and we can distinguish three 
classes of direct photon contributions differing by the overall heavy quark charge factor CF: 

Class A: CF = Eq. Comes from diagrams, like those in Fig. 3e-g, in which both photons 
couple to heavy QQ pairs. Despite the presence of mass singularities in contributions of 
individual diagrams coming from gluons and light quarks in the final state and from loops, 
the KLN theorem implies that at each order of the sum of all contributions of this class 
to cJdir is finite. Note that the first as well as the second terms in (3.1) arc also proportional 
to Bq and it is therefore this class of direct photon contributions that is needed for the 
calculation of Udir to be performed in a well-defined RS. 

Class B: CF = Cq. Comes from diagrams, like that in Fig. 3h, in which one of the photons 
couples to a heavy QQ and the other to a light qq pair. For massless light quarks this 
diagram has initial state mass singularity, which is removed by introducing the concept of 
the light quark (and gluon) distribution functions of the photon. The factorization scale 
dependence of the contribution of this diagram is then related to that of single resolved 
photon diagrams in Fig. 4a,c. 

Class C: CF = 1. Comes from diagrams in which both photons couple to light qq pairs, 
as those in Fig. 31. In this case the analogous subtraction procedure relates it to the 
single resolved photon contribution of the diagram in Fig. 4f and double resolved photon 
contribution of the diagram in Fig. 4h. The classes B and C are thus needed to guarantee 
the factorization scale (and scheme) invariance of the single and double resolved photon 
contributions to order a^a^. 

Because of different charge factors CF, the classes A, B and C do not mix under renormalization 
of as and factorization of mass singularities. As the diagrams in Fig. 3e and 31 give the same 
final state qqQQ, we should consider their interference term as well, but it turns out that it does 
not contribute to the total cross section £7(77 QQ)- 

3.3 Resolved photon contribution 




Figure 4: Examples of resolved photon diagrams involving the pointlike parts of PDF of the photon and 
the related direct photon diagrams. 
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The classes B and C of direct photon contributions of the order Q?c?g are indispensable to 
render the sum of direct and resolved photon contributions factorization scale invariant up to 
order a^a^. To see this in detail, let us write the sum of first two terms in (3.2-3.3) explicitly 
in terms of PDF and parton level cross sections 

a(i2)(M,/x) = 2as{li) jdxG{x,M) [^^^^(a;) + a,(^)c7^g(x, M, + 
•' i 

(^) JjdxdyJ2Qiix,M)qi{y,M) [ag^(ary) + a,(/x)(7g^(x2/, M, /x)] + 

alifi) J J dxdyG{x,M)G{y,M) [ag^(xy) + a,(/x)ag^(xy, M, /x)] + 

2al{fi) JJ dxdy^{x,M)G{y,M)af^{xy,M) (3.7) 

where qi runs over n f quark flavors and the factors of two and four reflect the identity 
of beam particles and equality of contributions from quarks and antiquarks inside the beam 
photons. Recalhng the general form of the derivative dcTres/dlnM^ 



dcJr, 



dlnM2 



dxWo{x,M) + dx 



dxdy 



J2 nix, M)Wg, (x, M) + G{x, M)Wg{x, M) 



G{x, M)Giy, M)WGGixy, M)+^ qi{x, M^y, M)Wgg{xy, M)+ 



+ 



T.{x,M)G{y,M)W^G{xy,M) 
using (2.4-2.6) and denoting = Q;s(/x), / = dJ/dlnM^ we find 

,2 



Wo{x,M) 



? ^kG\x)a^kx) + ^kf\x)Y,e-a^l{x,M) + 



(3.8) 



(3.9) 



W,,{x,M) = ^|47ra(J(x) + jdz[p^^l{z)a^l{xz) + ?,e',akf\z)a^^^ {xz)\ |+ • • • (3.10) 
Wg{x,M) = ^l^2T,a%:)+jdzPP^{z)a^^,{xz)^+--- (3.11) 

Wgg{x.M) = ^ |7ragi(x) +|dzp('g(z)agi(xz)} + • • • 
W,-,ix,M) = ^l^27ra^^{x) + 2ldzPS\z)al^{xz)^+--- 

W,g{x,M) = ^ |27rafJ(x) +/dzfp(°)(z)4)(x^) + P(°)(z)ag),(xz) 



(3.12) 
(3.13) 
(3.14) 



Only the lowest order terms on the r.h.s. of (3.9-3.14) are written out explicitly. All integrals in 
(3.10-3.14) go formally from to 1, but threshold behaviour of cross sections aij{xz) restricts 
the region to xz > Amq/W'^. 
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The factorization scale invariance of (3.7) requires that its variation with respect to InM^ 
is of higher order in than the approximation itself. There is no dispute that direct photon 
contributions of classes B and C are needed to guarantee this property. The question is which 
terms on the r.h.s. of (3.8) must vanish if the approximation is defined by (3.7). 

In the conventional approach both q{M) and G{M) are claimed to be of order a/ag and 
the approximation (3.7) thus of the order a^ag, implying that only terms up to this order must 
vanish in (3.8). This in turn means that the functions (3.10-3.11) must vanish to order and 
(3.12-3.14) to order Og respectively, which, indeed, they do ^. The fact that the expression on 
the r.h.s. of (3.9) does not vanish is of no concern in this approach as it is manifestly of the 
order aa^ and thus supposedly of higher order than (3.7) itself. 

If, on the other hand, we take into account that quark and gluon distribution functions 
of the photon behave as q{M),G{M) oc a, we see that Wq is of the same order a^a^ as the 
products qiWq., GWg and other integrands on the r.h.s. of (3.8) and must therefore also vanish 
for theoretical consistency of the approximation (3.7). This, in turn, necessitates the inclusion of 
class B direct photon contributions of the order aa^, like those in Fig. 4b,g, which provide the 
M-dependent terms the derivative of which cancels the first term in (3.8) involving the integral 
over Wq. Note that Wq^ in (3.10) receives contributions from the derivatives of both single and 
double resolved photon diagrams, proportional to cr^G" and figg, respectively. This fact reflects 
the mixing of single and double resolved photon contributions, which starts at the order a^al 
and is due to the presence of the inhomogeneous splitting terms in the evolution equations (2.4- 
2.6). For theoretical consistency of the sum of direct and resolved photon contribution up to 
the order aa^ only the lowest order double resolved photon contribution needs to be included. 

4. hh production at LEP2 

We now turn to the phenomenological analysis of bb production at LEP2, where the incoming 
leptons act as sources of transverse and longitudinal virtual photons, described by the fluxes 



PtM') = ^ ( - ^ ) , (4-1) 



n + {l-yf) 1 _ 2v4y 
27, \ y g2 Q. 

where stands for photon virtuality. Although the kinematic region of the LEP data includes 
photon virtualities up to moderate Q^, the cross section of the inclusive process 

e"'"e~ — s- e"'"e~65 -I- anything (4.3) 

is dominated by the production of the bb pair in the collision of two quasireal photons with 

very small Q^, typically (Q^) — 0.01 GeV^. For such small the cross sections of hard 
processes involving longitudinal virtual photons, which are proportional to Q^, arc expected 
to be negligible compared to those of transverse virtual photons. When talking about the 
production of bb in e+e~ collisions we shall always mean in association with the e+e~ pair, but 



^The latter condition is actually the same as for QQ production in hadron-hadron collisions. 
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for brevity of notation shall drop this latter specification, writing atot{^^^ bb). instead of 
o-tot(e"^e" e+e^bb). 

Although the data are available only for cross sections integrated over the whole phase space, 
we shall discuss the contributions dak{e'^e~ — ^ bb)/dW of individual processes as functions of 
77 collision energy W. The shapes of these contributions can altenatively be characterized by 
the functions 



i-W 

Fk{W)= / dw 

J2mb 



dukio^y- — hb) 
dw 



Gk{W) = 



dw 



bb) 



w 



dw 



(4.4) 



which quantify how much of a given contribution is located in the region up to a given W 
(Ffe(VF)) or above it (Gfc(l^)). As the available data are not copious enough to measure the 
differential distribution d(j{'s&e~ bb)/dW the theoretical analysis of the distributions (4.4) 
might allow us to invent a strategy how to separate the kinematic region of accessible W into two 
parts, each dominated by a particular contribution. The relative importance of the individual 
contributions as a function of W is determined by the ratia 



rk{W) = 



d(Tjk(e"' 



bb) I d(Ttot(e+e bb) 



dW 



dW 



(4.5) 



4.1 QED contribution 

The pure QED contribution to (7tot(e"'"e~ 

d(TQED(e+e- ^ 66) _ Ga^e^ AiW) 



bb) is given as 



dW 



ttS W 



1 + 



Ami 



8m^ 



where 



A{W) 



dydz6 



s 



yz 



i + (i-y)' 



1 + 



In 



^2^ 



1 + 
1-/3 



1 + 



4m^ 



2\ n 



(4.6) 



In 



^ max \ 



2 

max 



(1 



(4.7) 

results from convolution of of photon fluxes (4.1), integrated over the virtualities up to Q^ax- 
The convolution (4.7) can easily be performed analytically and the result inserted into (4.6). In 
Fig. 5 we display by the solid curve the result of evaluating (4.6) for nib = 4.75 GeV, ^/S = 200 
GeV and Q'^ax — ^ GeV. The distribution vanishes at the threshold W = 2m(, due to the 
threshold behaviour of the cross section (3.6), peaks at about W = 12 GeV and then drops 
rapidly off due to the fast decrease of both the photon flux (4.1) and (3.6). Integrating the 
distributions in Fig. 5 yields the values in the fourth column of Table 1. 

4.2 Leading order QCD corrections 

QCD corrections to pure QED expression (3.6) are of three types: direct (dir), single resolved 
(sr) and double resolved (dr). We shall first discuss the lowest order contributions to all three 
types of QCD corrections. As in the case of pure QED contribution, these corrections are given 
as convolutions of the photon flux (4.1) with the appropriate partonic cross sections. In all 
calculations u,d, s and c quarks were considered as intrinsic in the photon and Uf = 4 was taken 
in the expression for as(/i). 



- 10 - 




Figure 5: Left: The distributions dak/dW corresponding to the pure QED contribution together with 
three lowest order QCD contributions: single resolved (SR), double resolved (DR), and direct (DIR). The 
sum of all four contributions is shown as the upper dotted curve. Right: the comparison of pure QED 
contribution with the sum of three lowest order QCD ones. All curves were obtained for rrib = 4.75 GeV, 
^/S = 200 GeV and Q"^^^ = 4 GeV^ using GRV LO PDF of the photon and setting A^") = 0.27 GeV. 



Parameters 


QED 


LO QCD 


Total 


rub A(^) PDF 




DIR 


SR 


DR 




4.75 0.27 GRV LO 


1.27 


0.473 


1.415 


0.121 


3.28 


4.5 0.27 GRV LO 


1.40 


0.478 


1.746 


0.146 


3.77 


4.75 0.35 GRV LO 


1.27 


0.520 


1.542 


0.141 


3.47 


4.75 0.27 SASID 


1.27 


0.473 


0.904 


0.077 


2.73 



Table 1: The integrated cross sections cr(e+e~ bb,S) for y/S = 200 GeV and Ql^^x = 4 GeV^, 
corresponding to the distributions in Fig. 5. The renormalization and factorization scales fi and M we 
identified and set equal to uib. LO form of as{p) was used. All cross sections are in picobarns. 



4.2.1 Direct photon contribution 

The W dependence of the leading order QCD correction is given as the product 

of the convolution ^(M^) of photon fluxes and the lowest order QCD contribution as(^)o"j|](VK). 
At this order the direct photon contribution ct^|^ (//) , which comes form real or virtual emission 
of one gluon, is exclusively of class A. The function crj|^(PF/mft) has been calculated in, for 
instance, [22]. As it is just the first term in the series in positive powers of as(/u), the value 
of the renormalization scale /i in the argument of as(/x) is completely arbitrary. The resulting 
T^-dependence, evaluated for /i = rrib and shown in Fig. 5, is peaked even more sharply at 
small W than the pure QED contribution (4.6). This reflects the fact that the cross section 
c7j|](H^/mb) does not vanish at the threshold W = 2mb as does cr'^^iW/rri},). 
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4.2.2 Resolved photon contribution 

The leading-order single and double resolved photon contributions, were computed with HER- 
WIG Monte Carlo event generator, which implements the appropriate LO cross sections of the 
processes 

7 + G^6 + 6, (4.9) 
G+G^b + b, q + q^b + b, (4.10) 

where q = u,d,s,c stand for intrinsic quarks in the photon, and convolutes them with photon 
fluxes and PDF of the quasireal photon(s). In HERWIG the renormalization and factorization 
scales n and M are identified and set equal to an expression which is approximately equal the 
transverse mass /i = M = Mt = \J + . In LEP2 energy range the mean (Mt) depends 
weakly on W with, approximately, (My) ~ 7 GeV. 

Results of the calculations in which the LO GRV PDF of the photon, the LO expression for 
asifi) with A(^) = 0.27 GeV and rrih = 4.75 GeV were used, are shown in Fig. 5. As expected, 
the corresponding distributions are much broader than those of pure QED or LO QCD direct 
contributions. 

4.3 Comparison of individual contributions 

The comparison of the distributions da^/dW, Ffc(VF) and G^iW), corresponding to four individ- 
ual contributions, displayed in Figs. 5 and 6 and summarized in Table 1, reveals large difference 
in their shapes and magnitude. Specifically we conclude that 




Figure 6: Left: solid (dashed) curves show the partially integrated cross sections Fk{W) {Gk{W)) 
defined in (4.4) for QED and three LO QCD contributions. Right: The relative contributions rk{W) 
defined in (4.5) for the same four contributions. 



• The pure QED as well as the LO direct photon contributions peak at very small W and 
are basically negligible above ~ 50 GeV. For instance, the left plot of Fig. 6 shows that 
95% of the QED contribution comes from the region < 30 GeV. 
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• The onset of single as well as double resolved photon contributions is much slower, but 
these distributions are, on the other hand, markedly broader. 



• The double resolved photon contribution is negligible everywhere. 

• The pure QED and single resolved photon contributions are of comparable size and to- 
gether make up about 85% of the total integrated cross section, 

• up to about 2± 30 GeV, datot/dW is dominated by pure QED contribution, whereas for 

> 30 GeV, QCD contributions take over. 

The numbers given in Table 1 correspond to standard colored quarks with fractional electric 
charges. In [16] the excess of data over standard theoretical calculations is interpreted as evidence 
for Hahn-Nambu integer quark charges. Applied to the case of b-quark, the author of [16] argues 
that the correct way of calculating the charge factor in (4.6) is not the usual 3e^ = 1/27, but 



where the sum runs over the three Hahn-Nambu integer b-quark charges , which are 0, 0, —1 
respectively. The results is thus 9 times more than in the standard calculation. I think his 
argument for first summing in the formula (4.11) over the quark colours and then taking the 
fourth power is wrong ^, but I mention it here because were the author of [16] right, the whole 
discrepancy would come from the region of small W, where QED contribution dominates. 

On the other hand, were the light gluino production [15] responsible for the observed excess, 
the latter would have to come from the region of W dominated by the double resolved photon 
contribution. Although the energy dependence of the gluon-gluon fusion to gluino-antigluino 
may be slightly different than those of G + G QQ ot qq ^ QQ, it is clear that the basic shape 
of the T^-distribution is given by the convolution of the photon fluxes (4.1-4.2) and the gluon 
distribution function of the photon, which are the same in both types of processes. 

The above observations underline the fact that in order to pin down the possible origins of 
the excess of the integrated cross section atoti^~^^~ e^e~bb) over theoretical calculations, it 
would be very helpful if the data could be separated at least into two subsamples according to 
their hadronic energy W, say W <30 GeV and W >30 GeV. 

4.4 Dependence on input parameters 

The magnitude of the contributions discussed in the preceding subsection depend, beside the 
e"'"e~ cms energy -v/S", on a number of input parameters: the numerical values of mi,, Aq^j^, Q^ax^ 
the selection of PDF and the choice of the renormalization and factorization scales and M. In 
all the calculation reported above we set ji = M = mi,. The central calculation was performed 
for y/S = 200 GeV, Q^„^ = 4 GeV^, = 4.75 GeV, A^'') = 0.27 GeV using the GRV LO PDF 
of the photon. To see the sensitivity of the LO results to these assumptions we varied some of 
these parameters: 

^The correct value of the charge factor in (4.6) for the Hahn-Nambu integer charge b-quark equals 1 and would 
thus yield even bigger enhancement than that suggested in [16]. 




(4.11) 
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• mi, was lowered to m;, = 4.5 GeV, 
was increased to 0.35 GeV, 



• GRV set of PDF of the photon was replaced with that of Schuler-Sjostrand set SASID. 

The choice of Q^ax — 4 GcV^ corresponds roughly to the usual cuts imposed on the LEP2 data 
and could therefore be also adjusted to specific conditions of a given experiment. 

The results of the calculations of crtot(e''"e~ — e+e~66), corresponding to different sets of 
input parameters specified above, are listed in Table 1. Lowering rrib increases all four contribu- 
tions, as docs, except for the pure QED one, increasing A^''). SASID PDF yield markedly lower 
results for single and double resolved photon contributions. It is clear that varying the input 
parameters within reasonable bounds docs not bring the sum of lowest order QED and QCD 
calculations significantly closer to the data. 

5. Can the next-to-leading order QCD corrections solve the puzzle? 

With the sum of lowest order QED and QCD contributions to cj^j^(e^e^ e^e^bh) way below 
the data one might expect that the NLO contributions could at least partly bridge the gap 
between data a theory. 

5.1 Direct photon contribution 

The the sum of the second and third terms in (3.1) can be written, suppressing the dependence 
of cr^j|. on the ratio W/mi,, as 



dir 



f^dir^^lM) [1 + ri{ii/mb)as{n)] . (5.1) 



Note that W^^cr^|^ is a unique function of the ratio W/mb and the NLO coefficient ri(W,mh, fi) 
can be written as a function of W/nib and nih/lJ- The first term in (5.1) is a monotonous 
function of the renormalization scale fi, spanning the whole interval between zero and infinity. 
As emphasized in Section 3.2, one needs to include at least the term OjC^j^ to make the expression 
(5.1) of genuine next-to-leading order. Only the class A of order a^a^ direct photon contributions 
is needed for this purpose. The renormalization scale invariance of cr^^^ implies the following 
general form of ri: 

n{W/mb,fx/mb,RS) = |° In-^ - p{W/mb), (5.2) 

where p{W/mb) is a renormalization scale and scheme invariant [17], which, however, depends 
beside the ratio W/mi, also on the numerical value of the ratio mi)/k^^\ It can be evaluated 
using the results of a calculation in any renormalization scheme RS 



p(M/^/mb,mb/AW) = ^ In ^ - ri(W^M, 1, RS) (5.3) 



and its numerical value governs basic features of the scale dependence of (5.1): 
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Figure 7: a) The renormalization scale dependence of the leading (solid curve) and next-to-leading 

order contributions to the generic quantity (5.6) for different values of p. The dashed curves correspond 
to p > 0, the dotted ones to p < 0, the dash-dotted to p = 0; b) the same as in a) but plotted in a linear 
scale of /Lt; c) graphical representation of (5.6) in three different renormalazation schemes and for p > 0; 
d) the shape of the NLO expression (5.6) as a function of ag. 



p > 0: the NLO approximation (5.1) considered as a function of /x exhibits a local maximum, 
where da^^^ /d/j, = and where the prediction is thus most stable. This point, preferred 
by the Principle of Minimal Sensitivity [17], is also very close to the point for which ri = 0, 
which is selected by method of Effective Charges [21]. The value of cr^^'^ at this point 
is proportional to 1/p implying very large NLO corrections for small p. Inserting the 
appropriate numbers for n/ = 4, rrib = 4.75 GeV and = 0.27 GeV, we get 

p{W/mb) = 3.88 - ri{W/mb, 1,MS). (5.4) 

The NLO coefficient ri{W/mb, 1,MS) does not therefore have to be outrageously large to 
get small, or even negative p\ 

p < 0: cr^ir*^ becomes a monotonous function of p, similarly to cr^^. In fact, taking the 
derivative of o'^^'~' with respect to p one finds that for p < it is actually even steeper 
than cT^j^, given by the first term in (5.1). Consequently, for negative p going to the NLO 
does not improve the stability of the calculation, but quite on the contrary! 
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The above features are straightforward to see setting Pi = 0. This technical assumption simphfies 
the relevant formulae significantly, but nothing essential depends on it. Setting c = in (2.2) 
allows us to write Q;^j(/i) explicitly as 

47r 

which, inserting this expression into (5.1) and taking into account (5.2), gives 

'^dir° = ^S«.(/^)[2-p«.(/x)], (5.6) 

In Fig. 7 we plot the dependence of the generic NLO quantity (5.6) for several values, positive 
as well as negative, of p and in different renormalization schemes. Several conclusions can be 
drawn from this figure: ^ 

• For p < <7^ir*^, displayed in Fig. 7a, is a steeper function of fi than a^P because the 
second term in (5.1), equal to 1 — pas{p.) is monotonously decreasing function of p. For 
positive p, on the other hand, (5.6) exhibits a local maximum at a^"^ = l/p- 

• For negative as well as positive values of p, cr^^'^ ocl/ln|Uas;u— >-oo. For negative p 
this implies that there is no region of local stability except for the trivial one at /x = oo. 
However, when plotted on a linear scale of ^ in a limited interval the weak logarithmic 
does, as illustrated in Fig. 7b, fake the local quasistability. 

• The curve representing (t^^'^ depends, as shown in Fig. 7c for p > 0, on the chosen 

renormalization scheme. Setting fi equal to some "natural" physical scale Q therefore 
does not resolve the renormalization scale ambiguity as in different RS we get different 
results. However, although also the position of the local maximum depends on the choice 
of the RS, the value of the NLO approximation (5.6) does not! The same holds for the 
intersection of the LO and NLO curves. 

• Instead of varying both the renormalization scale and scheme, which is legitimate but 
redundant, we may use the couplant ag itself for labeling the different predictions of (5.6). 
Instead of an infinite set of curves describing the //-dependence of a^^^{p,, RS) in different 
RS, we get for each p a single curve displayed in Fig. 7d. 

We wish to emphasize that whereas there are natural physical scales in hard collisions, there is 
nothing like the "natural" renormalization scheme. There are only two general lines of arguments 
for choosing the renormalization scheme: either one looks for the maximum local stability or 
smallest order corrections. In the first case one is led to the Principle of Maximal Stability [17], 
in the second to the method of Effective Charges [21]. As illustrated by the situation for negative 
p, there may, however, he cases when none of them can be applied. 

As the term in (5.1) has not yet been calculated, we cannot associate the class A 
direct photon contribution, given by the first term in (5.1), to a well-defined renormalization 



^These conclusions are well know in the context of perturbative quantities of the form (5.1) depending on 

the renormalization scale only. We recall them because they will be used in the analysis of the resolved photon 
contribution, where the interplay between the renormalization and factorization scales complicates the situation. 
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scheme. Moreover, as the magnitude of the NLO corrections in (5.1) is determined by the ratio 
a'^^{W/mb, 1) / al^^ {W / rrib) of two functions of W/rub, which may depend on W/nib in different 
ways, the coefficient ri may be very large even when both the numerator and denominator are 
on average of comparable and small magnitude. The size of NLO correction may also depend 
sensitively on the ratio W/nib- All this indicates that without the knowledge of class A direct 
photon contribution of the order a^a^ , we cannot make a meaningful estimate of the importance 
of higher order corrections to (5.1). 

5.2 Resolved photon contribution 

As in LEP2 energy range the double resolved photon one is numerically negligible, only the single 
resolved photon contribution will discussed in detail below. As shown in Fig. 5, the spectrum of 
the contributions dagr/dW peaks at about W = 30 GeV with the mean value (W) = 65 GeV. 
The properties of the measures cross section (TjjJ:'*-'(e^e~ — > e'^e~bb) will therefore be determined 
primarily by those of o'^^^i'yj —>■ bb) in the energy range 30 < < 65 GeV. All the results 
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Figure 8: a) The scale dependence of the conventional NLO approximation o'^'"*-'(VL, M, M) for = 40 
GeV (solid curve) together with the contributions of the "fG (dashed) and ^q/q (dotted) channels. The 
LO approximation is shown for comparison by the dash-dotted curve, b) (j^^'^{W, M, kM) (solid curve) 
and crLO(VL, M, kM) (dash-dotted) for VL = 40 GeV and three values of k = 0.5, 1.2. 

shown below for the single resolved photon contribution 



(J, 



NLO/ 



{W,M,fi) = 2a,(/x) J dxG{x,M) ^a^'^(x) + a,(/x)a;^^(x, M, /i) 
Aalifi) [ dxY,Q^i^,M)a!^^^{x,M) 



.(2), 



+ 



(5.7) 



are based on the formulae for the partonic cross sections a^j^ as given in [23]. Even if the 
reader does not agree with our claim that the approximations employed in [7-9] do not consti- 
tute complete NLO approximation, it is certainly important to understand quantitatively their 

(2) 

renormalization and factorization scale dependence. Because the expressions for a^^ as given 
in [23] correspond to fi = M, we have restored its separate dependence on jj, and M by adding 
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to o"^,q(.t, M, M) the term (/?o/47r)(j^^ ln(/i^/Af^). Note that for each vahie of factorization scale 
M the expression (5.7) has, as far as the //-dependence is concerned, the form of the NLO ex- 
pression (5.1). In all calculations the GRV HO set of PDF of the photon and A^^^ = 0.27 GeV 
were used. 

We first follow the conventional procedure and set M = fi. The resulting (common) scale 
dependence of the expression (5.7), together with those of the quark and gluon contributions to 
it, are shown in Fig. 8a for = 40 GeV. Overlaid for comparison is also the LO approximation, 
given by the first term in (5.7). We note the different scale dependence of the jG and 75 channels, 
the latter turning negative for M > 6 GeV, but the most important observation concerns the 
fact that the conventional NLO approximation (5.7) is a monotonously decreasing function of 
the scale. Moreover, it falls off even more steeply than the LO expression! In other words in 
going from the leading to the next-to- leading order the sensitivity to the scale variation increases, 
rather than decreases, as one might expect (and hope)! Recalling the discussion in Section 




20 



Figure 9: The surface and contour plots of a.^^ (W^- /") for W = 40 GeV. 



5.1 one should, however, not be surprised. To check how much this feature depends on setting 
exactly 11 = M, we plot in Fig. 8b the scale dependence of cj^^*-'(H^, M, /x = kM) for standard 
choices of At = 0.5, 1, 2. Clearly, the above conclusion is independent of k in this range. 

The steep and monotonous scale dependence of cr^^^(VF, M, = kM) is a clear warning 
that the conventional NLO approximation is highly unstable. To see what happens if we relax 
the usual but arbitrary identification = kM we plot in Fig. 9 the surface and contour plots 
representing the full M and ^ dependence of cj^'^^(VF, M, as given in eq. (5.7). Contrary to 
analogous process in antiproton-proton collisions [12], it does not exhibit a saddle point, where 
the derivatives with respect to both M and /i vanish, but Fig. 9 seems to indicate some sort of 
stability region at large scales, say for M > 10 GeV, /x > 20 GeV. This impression is, however, 
misleading as becomes clear if we plot in Fig. 10 the slices of the surface plot in Fig. 9a along 
both axis and recall the discussion of Section 5.1. For each fixed value of the factorization scale 
M the expression (5.7) has a form of the NLO expression as far as the renormalization scale 
/i is concerned. Comparing the curves in Fig. 10a with those of 3a we see that for M < 4.2 
GeV a^^^{M,ii) corresponds to negative p in (5.6) and thus exhibits no local stability point. 
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Figure 10: The renormalization scale dependence of afj}"'~'{W, M, /i) for fixed values of the factorization 
scale (a) and vice versa: the factorization scale dependence of af^^''^ {W, M , ^) for the same set of fixed 
values of the renormalization scale fi (b). All calculations correspond to = 40 GeV. In a) the ordering 
from above of the curves at ^ — 30 GeV corresponds to M — 30,16,10,7,4.75,3,2,1 GeV, in b) the 
curves correspond at M = 0.7 GeV to the same sequence from below. 



For higher M the local maximum in the /i-dependence of o"^ (M, /i) exists at the associated 
^J-maxiM). The M-dependence of a^j^'^{M, fj,rnax{M)), shown in Fig. 10a by the dotted curve, 
is, however, even steeper that those at fixed M. The above plots and conclusions concerned 
the results at one typical value of W, but their essence holds for the whole interval relevant for 
LEP2 data. 

We thus conclude that in the energy range relevant for LEP2 data the renormalization and 
factorization scale dependence of the conventional NLO calculations of single resolved photon 
contribution to the total cross section cJtot(77 — > bb) exhibits no stability region, either as a 
function of the common scale fi = kM or as fully two dimensional function of /U and M. 

6. Summary and Conclusions 

We have argued that in order to understand the origins of the excess of LEP2 data on bb 
production in 77 collisions over the current theoretical calculations, two ingredients are needed. 

On the experimental side, the separation of data into at least two bins of the hadronic 
energy W, say W < 30 GeV and W > 30 GeV, could be instrumental in pinning down the 
possible mechanisms or phenomena responsible for the observed excess. 

On the theoretical side, the evaluation of the direct photon contribution of the order a^a^ 
is needed to provide the terms that would make the existing theoretical expressions of genuine 
next-to-leading order in a^. In their absence, the existing NLO calculations are highly sensitive 
to the variation of renormalization and factorization scale and thus inherently unreliable. 
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